Quantifying entanglement with covariance matrices 
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Covariance matrices are a useful tool to investigate correlations and entanglement in quantum 
systems. They are widely used in continuous variable systems, but recently also for finite dimensional 
systems powerful entanglement criteria in terms of covariance matrices have been derived. We show 
how these results can be used for the quantification of entanglement in bipartite systems. To that 
aim we introduce an entanglement parameter that quantifies the violation of the covariance matrix 
criterion and can be used to give a lower bounds on the concurrence. These lower bounds are easily 
computable and give entanglement estimates for many weakly entangled states. 

PACS numbers: 03.67.-a, 03.65.Ud 



I. INTRODUCTION 

Entanglement is a central resource in quantum infor- 
mation processing and many works are devoted to its 
characterization [l|, 0, [![• One line of research is the 
derivation of entanglement criteria which should detect 
the entanglement also of weakly entangled states. A dif- 
ferent line of research tries to quantify the entanglement 
via so-called entanglement measures As most en- 
tanglement measures arc defined via complex optimiza- 
tion procedures, they are often difficult to compute and 
therefore one tries to give at least lower bounds on them 

an. 

Covariance matrices (CMs) of local observables are a 
widely used tool to study correlations in continuous vari- 
able systems, such as coupled harmonic oscillators or 
modes of light Q. Moreover, many entanglement cri- 
teria for these systems are formulated as conditions on 
CMs [H . Recently, it has been shown that CMs are also 
a useful tool for the investigation of entanglement in dis- 
crete systems, such as polarized photons or trapped ions 
Hi OH EH- Indeed, in Refs. PUEI a so-called covariance 
matrix criterion (CMC) has been established, which al- 
lows to detect many weakly entangled states, which are 
not detected by other criteria. 

In this paper we show that CMs can be used not only to 
detect but also to quantify entanglement in discrete com- 
posite quantum systems. To that aim, we define an en- 
tanglement parameter £(g), that quantifies the violation 
of the CMC for discrete systems. Our construction is in- 
spired by a similar definition of an entanglement param- 
eter for continuous variable systems in Ref. [l2|]. While 
it remains unclear to which extent £{g) is a direct entan- 
glement measure, we will see, however, that £(g) gives a 
lower bound on the concurrence, which is a widely used 
entanglement monotone. 

In detail, our paper is organized as follows: In Sec- 
tion II we recall the CMC and define our entanglement 
parameter £. We also demonstrate that existing results 
about the CMC give directly lower bounds on£. In Sec- 
tion III we consider general properties of the the param- 



eter £ . We show that is is convex and invariant under 
local rotations, but we present an example, where £ in- 
creases on average under LOCO The physical reason be- 
hind this example is the fact that there are entangled 
states, which cannot be detected by the CMC, however, 
they are detected after suitable local filtering operations. 
In Section IV we investigate, how £ can be computed 
for different types of states. We give explicit formulas 
for pure states, and also show how to compute £ for the 
family of mixed, Schmidt-correlated states. In Section 
V we show how £ can be used as a lower bound on the 
concurrence. This delivers non-trivial bounds on the con- 
currence for many weakly entangled states. Finally, we 
conclude the paper and give technical calculations for 
some of our theorems in the Appendix. 



II. DEFINITION OF THE ENTANGLEMENT 
PARAMETER 

In this section we introduce a function based on the 
CMC that can be used to estimate the amount of entan- 
glement in a given quantum state. To start, let us fix 
the notation and introduce the quantities, which we are 
going to work with. We consider quantum states g over 
a finite-dimensional bipartite Hilbert space Ha ®Hb, 
where cIa = dim(7i J 4) and ds = ctim(Hg) denote dimen- 
sions of corresponding local spaces. Physical observables 
are described by Hermitian operators. For our purpose 
we choose a complete set of orthogonal observables Ai on 
Ha with i = l,...,d 2 A and Tr (Ai A j) = Sij and a similar 
set Bj for Hb- We will refer to them as local orthogonal 
observables, an example for dA = 2 are the (appropri- 
ately normalized) Pauli matrices and the identity. Then, 
we can consider observables on Ha ®>Hb defined by 



{M a } = {A, 



3=d\ 



1,1® Bj}, i = l, 
1,- 



'A- 



,4 



4, 



(1) 



which then also obey Tx{M a Mp) = 6 a p. 

The main object of our studies will be covariance ma- 
trices (CMs). A CM of a given bipartite state g is defined 
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by the following entries 

l(Q)*[) = \{M a M + MpM a ) e - {M a ) e {M ) e . (2) 

Choosing the observables as in Eq. |T]) one can write the 
CM in a handy block form 

7={c t b)' (3) 

where A = "/(qa, {Ai}), B = ^(qb, {-£>»}) are CMs of re- 
duced density matrices and Cy = (Ai®Bj) e ~(Ai) g (Bj) e 
denote correlations between the two parties. 

Before introducing the function that we are going to 
use for entanglement quantification let us state the co- 
variance matrix criterion (CMC). For that, recall that a 
state is separable, if it can be written as a convex combi- 
nation of product states, i.e. g = J2 k Pk\ak){ak\ © \bk){bk\ 
with some probabilities pk- Then we have: 

Theorem 1 (Covariance matrix criterion). Let g be a 

separable bipartite state. Then there exist pure states 
IV'fcMV'fcl i> n Ha and |<^fc)(<^fc| in Hb and convex weights 
Pk such that if we define ka = '^2kPkl{\4'k}{'4'k\) and 
k b = J2kPkl(\4>k)((f>k\) the inequality 

i(e,{Mi}) >k a ®k b ^ > fa 1 (4) 

holds. This means that the difference between left and 
right hand side must be positive-semidefinite. If there 
are no such ka,b then the state g must be entangled. 

The proof of this statement can be found in Ref . [l(| • 
The main task for applying the CMC is its evaluation, 
that is, the characterization of the matrices ka and kb- 
For this, several corollaries of the CMC have been derived 
in Refs. [13, [HI- As we will use them later, we present 
some of them here, but without any proof. For simplicity, 
we only consider the case d = dA = ds- 

Proposition 2 (CMC evaluated from traces). Let g be 
a state with CM 7 as in Eq. (0). Then if g is separable, 
we have 

d 2 d 2 
22K|C|)< (j2 Aii - d + 1 ) + (Y, B "- d+1 ) 

i=l i=l 
= [l-Tr(g A )] + [l-Tr(gl)}, (5) 

If this inequality is violated, then g must be entangled. 

Proposition 3 (CMC and the trace norm of C). Let g be 

a state with CM 7 as in Eq. (Q). Then if g is separable, 
we have for the trace norm of C 

\\C\\l<[l-Tr{gl)}[l-TriQ%)}, (6) 

If this inequality is violated, then g must be entangled. 



In order to define our entanglement parameter £, let us 
reformulate the CMC in a slightly different way. Imagine 
some state g is detected as entangled by the CMC. On 
the one hand there exist no ka and Kg as above such 
that 7(g) — ka © kb > 0. On the other hand we can find 
surely k\ and k b and some number t e € [0, 1] such that 
7(0) ~ t e K e A © K e B is again positive semidefinite: 

7(p e ) - t e K e A ®k b >0. (7) 

In the worst case, we can fulfill this inequality by choosing 
t e = 0. For a state that is not detected by the CMC (e.g. a 
separable state) the parameter t can be chosen to be at 
least one, or even larger than that. 

Implementing this idea in Theorem [T] results in an al- 
ternative formulation of the CMC: 

Theorem 4 (Parameterized CMC). Let g be a bipartite 
state. Assume that we choose pure states \ipk)(tpk\ on Ha 
and \(f> k )((f>k\ on H B such that k° a = J2kPkl{\i'k)(tpk\) 
and k° b = J2kPkli\4'k}(4 l k\) are optimal in the sense that 

7 - t K° A ®K° B >0, (8) 

for some < t a < 1, but 

7 — tKA © kb ^ 0, for all t > t a and all ka, kb- (9) 

Then if the state g is separable there exist k° a and k° b 
such that 

max{i < 1 : 7 - tn° A © K° B > 0} = 1, (10) 

otherwise the state is entangled. 

This leads to the idea, to use for entangled states the 
parameter t Q as an entanglement parameter. More pre- 
cisely, we can define: 

Definition 5 (Entanglement parameter). Let g be a bi- 
partite quantum state with CM 7(g). We define a func- 
tion V(g) as 

V(g)= max {t < 1 : 7(e) - tK A © k b > 0}. (11) 

t,KA,KB 

The entanglement parameter £{g) is then defined as 

S( Q ) = 1-V(q). (12) 

The parameter £(g) vanishes for separable states and 
is larger than zero for all states that are detected by the 
CMC. This function £ (g) is the main topic of study in 
this paper and, as we shall see later, can be used to quan- 
tify entanglement in quantum states. A similar function 
has been already used to quantify entanglement in infi- 
nite dimensional systems, namely Gaussian states 
there this parameter turned out to be an entanglement 
monotone for special operations on special states. 

Interestingly, using the parameterized version of the 
CMC (Theorem 2]) and Propositions and [3] one can 
immediately give a lower bound £{g). We can formulate: 
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Proposition 6 (Bounds on £{g)). Assuming that d 
d>A = d-B we have in the situation from above that 



Tr(gl)+Tr(g%) + 2Tr(\C\)-2 



2d -2 



(13) 



and 



{{Tr(g\)-Tr(g%W + \\C\\i}- (14) 



Proof. For the first calculation as in Ref. fll| 

gives a parameterized version of Proposition^ and results 
in 2Tr(|C|) < Tr(A + B - t(n A + k b )). Using Tr ( 7 (g)) = 
d - Tr(g 2 ) (see Ref. [U) gives 

2rf-Tr(g^)-Tr(g|)-2Tr(|C'|) 
" 2^~2 (15j 

and finally Eq. (fl3|) . Eq. (fl4|) can be also directly derived 
from Eq. ([6]) and from the calculations in Ref. Ill] . ■ 



III. PROPERTIES OF THE ENTANGLEMENT 
PARAMETER £ 

In this section we investigate general properties of the 
function £{g). Since the function £(g) should be used 
to quantify entanglement in a given quantum state, two 
of the properties that have to be fulfilled are that it is 
convex and does not change under local unitary transfor- 
mations. Indeed, this is the case: 

Lemma 7 (Convexity and invariance under local uni- 
tary transformations). The entanglement parameter £{g) 
is invariant under local unitary transformations and is 
convex in the state, that is for g = pg\ + (1 — p)Q2 we 
have that £{g) < p£{gi) + (1 — p)£(g 2 ). 

Proof: The invariance under local unitary transforma- 
tions follows simply from the fact that the CMC is invari- 
ant under such transformations [l(| [HI ■ In more detail, 
such transformations map a set of local orthogonal ob- 
servables to another set of local orthogonal observables, 
and the CMC does not depend on the choice of the ob- 
servables. 

Concerning convexity, it is sufficient to prove the con- 
cavity of V(g), i.e. that for any state g = pg\ + (1 — p)g 2 
the inequality V(g) = t > pi\ + (1 — p)i 2 = t' holds, 
where i\ = U(gi) and i 2 = V(g 2 ). 

To prove this we exploit the connection between 
the CMC and local uncertainty relations (LURs) [l3|. 
V(g) = t implies that the parameterized CMC crite- 
rion is fulfilled and there exist ka, hb and t such that 
■y(g) — tKA © kb > 0. According to the Proposition V.2 
in Ref. [ll| this means that if we take arbitrary local ob- 
servables on Alice's and Bob's side Ak © 1 and 1 <%> Bk 



such and define positive constants XJ a — min e 8 2 (Ak) 
and Ub = min e ^ fc 8 2 {Bj t ) then 

Y,S 2 (A k ® l + l®B k ) g >t(U A + U B ). (16) 



Therefore it suffices to show that t' fulfills the last in- 
equality as well. Due to the concavity of the variance we 
can write 

52 ( A k ® 1 + 1 ® B k ) g > pJ2 52 ( A k ® 1 + 1 ® B k ) ei 

k fc 

+ {\-p)Y J 5 2 {Ak®t + t®B k ) g2 . (17) 

fc 

Since the states Q\ and g 2 both fulfill the CMC with the 
parameters t\ and ti we can write 

pY, 6 2 {A k ® 1 + 1 ® B k ) ei + (1- P )J2 S 2 {A k ® 1+ 
fc fc 

+ l®B k ) e2 > [t lP + t 2 (l-p)] (Ua + Ub). (18) 

Note that t is defined as maximal value of all possible t. 
Using (JT7J) and (jT5J) this finishes the proof. ■ 

A further important property of entanglement mea- 
sures is they do not increase under local operations as- 
sisted with classical communication. This condition can 
be demanded in two different forms (see Refs. fll. 13. [l4| ) : 
Minimally, one requires that if g arises from g via some 
LOCC transformation, then E(g) > E(g) holds. Often, 
however, a stronger condition is required and fulfilled, 
namely that E(g) should not increase under LOCC oper- 
ations on average. This means that if an LOCC protocol 
maps g onto some states gi with probabilities pi, then 



E(g)>Y,P* E (Qi)> 



(19) 



should hold. 

In the following, we will show by an example that £{g) 
can increase on average under LOCC operations. This 
does not exclude a priori the usability of £{g) as an entan- 
glement monotone (since the minimal requirement might 
still hold), however, it is a hint that £(g) might not be 
an entanglement measure. As we will see later, however, 
£ (g) can be very useful to derive lower bounds on the 
concurrence for mixed states. 

Lemma 8 (Increasing one average under LOCC). There 
exists a two-qubit state g and an LOCC-protocol, such 
that £ increases on average from zero to a positive value 
under this protocol. 

Proof. We prove the statement by providing an explicit 
example of a two-qubit state, which can be found numer- 
ically. The idea to find such an example is as follows: We 
consider a family of states that was already intensively 
investigated in Refs. [ll], EH- Within this family one can 
find pairs of states g and g' with the same covariance 
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matrix but where g is entangled, while g' is not. Hence, 
g cannot be detected by the CMC criterion, and £(g) has 
to vanish. 

It was shown in Refs. pUdU, however, that after an 
appropriate filtering operation 



g >-> gmt = F A ® F b qF\ <g> F B 



(20) 



any entangled two-qubit state can be detected by the 
CMC. Hence £{gm.t) > and the filtering operation will 
give rise to the desired LOCC operation. 

To be more concrete, a numerical example of the afore- 
mentioned state g is 



/ 0.48508 0.02094 

0.33 

0.00067 

V 0.02094 0.18425 



(21) 



which is not detected by the CMC (see but which 
is clearly NPT and hence entangled. The corresponding 
filter operations are 



F A = 
F B = 



0.16457 

0.98637 

0.96526 

0.26128 



(22) 



The final state after gfm filtering will be 
F A ® F B g'F A ®F B 



gmt 



Tr (F A ®F B g'F A 
/ 0.47636 



F B ) 


0.02375 
0.02364 

V 0.03336 




(23) 



and is detected by the CMC, hence S(Qmt) > 0. Since g 
is not detected, we have £{g) = 0. 

Using the filter operations F A and F B we can now con- 
struct a POVM type of measurements for Alice and Bob. 
The complementary operations are given by 



F% = (l-F A F A )i 
F B = (l-F B F B )i 



0.98637 

0.16457 

0.26128 

0.96526 



(24) 



With this operations we establish LOCC protocol with 
four different outcomes 

Qx = gfiit with probability p\ = 0.02570, 
02 with probability pi = 0.17629, 
g3 with probability p^ = 0.46200, 
qa with probability p^ = 0.33601. (25) 

Important for us is the fact that applying this pro- 
tocol to a state with £{g) we achieve a state such 



that £(f?fiit) with non-zero probability. Therefore = 

< Pi£(ffi)i ancl £{q) increases on average un- 

der LOCC. ■ 

Note that for the provided example one can check the 
separability of the state g = ^ li Pigi as this state has 
a positive partial transpose and is therefore separable. 
Consequently, the protocol given is not a counterexample 
to the LOCC condition of the first kind. 



IV. EVALUATION OF £(g) FOR PURE AND 
SCHMIDT-CORRELATED STATES 

In this section we compute £{g) for pure states and a 
family of mixed states. We start with the case of two- 
qubits. Then, we generalize it to d-dimcnsional systems. 



A. Pure states of two qubits 

Using the relations that can be found in Appendix A, 
it is straightforward to calculate the CM of a two-qubit 
state \ip) = V%"|00) + with Ai + A 2 = 1. The 

CM will have the familiar block form 



( C T B 



(26) 



with 



A = B 



V 



\ 2 

1 - A l 
A1A2 






-A1A2 







c 



f\i-X\ 

-A1A2 



V 



— Ai A2 
A 2 - Ai 







VA1A2 




(27) 



VA1A2 / 



The next step in the calculation of the parameter t and 
therefore of the function £{g) is to find the optimal 
k a ® k b . In the two-qubit case we first guess the correct 
solution and the prove its optimality. 

To construct the matrix k a ® n B we take two prod- 
uct states |00)(00| and |11)(11| and get k a © n B = 
idiag{0, 0, 1, 1, 0, 0, 1, 1}. Then we calculate the V{\ip)) 
from the condition 7 — tn A © k b > 0. This matrix is 
positive iff 1 — t > 2y/X 1 \2 and therefore for any 



t<l- 2^/M\2 



(28) 



we can find k a and n B such that 7 — tn A © k b > holds. 

Note that taking some particular expansion for k A © 
k b , strictly speaking, does not provide any information 
about the entanglement, except for the case when we are 
able to find k a ® n B such that 7 — tn A © k b > for 
some t > 1. Then the state is not detected by the CMC 
and £ (g) = 0. However, we can use the Proposition [5] to 
prove the following: 
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Lemma 9. The upper bound on the parameter t for two 
qubits provided in Eq. 128\) is tight. 



Proof: Directly applying the relation (fT5j) to the two- 
qubit case we have t < 1 — 2^/\~[X^, which coincides with 
(|28|) and therefore gives an optimal bound on parameter 
t. Indeed, on the one hand, it follows immediately from 



(|28| that if t < 1 — 2v%A2 then we can find a decom- 
position ka® kb such that 7 — tKA © «b > holds. On 
the other hand, the condition (fTSjl implies that for all £, 
with i > 1 — 2VA1A2 and for all re^ and the relation 
7 — tKA © Kb ^ holds. ■ 

According to the last Lemma the function £(\ip)) can 
be calculated exactly for two-qubit pure states as 



B. Pure states of two qudits 



(29) 



To estimate the parameter t for a pure state of two 
d-level systems, we follow the same strategy as in the 
two-qubit case and take the states \kk) (kk\ for the decom- 
position of ka © Kb m order to derive the upper bound 
on the parameter t. We make the ansatz 



ka © Kb 



ti 

E 



Pil{\H)) 



(30) 



with some probabilities pi . 

The positive semi-definiteness of the matrix 7 — tKA © 
kb then implies the positive semi-definiteness of 2 x 2 
blocks of the type 



yy _ 

^2x2 — 



A,; + Xj - t(p.i + Pj) 



±2 V /A~A, 

for all i < j. Therefore, if for all i < j 



±2y%X~ 

A; + Xj - t(pi - 



Pj) 



t < 



Pi + Pj 



(31) 



(32) 



holds, then we can find ka and kb such that 7 — tKA © 
kb > holds. To achieve the goal and calculate the 
function £(\ip)) we need to prove that the choice of the 
expansion of the ka © Kb in Eq. (|30|) was optimal. 



Lemma 10 (Optimality of the decomposition). The op- 
timal expansion for ka® kb can always be written in a 
form of the Eq. 



^opt 



© K 



opt 



E 



Pil{\H))- 



(33) 



Proof. First, we show that for pure states in Schmidt 
decomposition 7(|V ; )) ~ t- K A © Kb > is equivalent to 
j(\tp)) — tK®K > 0, for some re, which can be found explic- 
itly. This re can be constructed by choosing the product 



states in a proper way. Indeed, note that since the CM 
of a state in Schmidt decomposition is symmetric with 
respect to the interchange of the parties (A <-» B) the 
relation 7(1^)) — tKB © ka > must hold as well. So let 

K A = J2k=iPkl(\ a k)(ak\) and k b = J2k=iPkl(\ b k)(h\)- 
Then we have 

7(IV>» 



(ka © k b + k b © Ka) > 0. 



(34) 



Since ka © kb + kb © ka = ka © ka + kb © kb, the 
appropriate choice of the product states is 



\Vk) 



. , K (for ka © re^), 
1,...,2K (for k b © Kb)- 
(35) 

Hence we have t( I V 7 ) ) — tre © re > 0, with 



\a t ) © \di), i = 1, 
\h) © \h), i = K 



2K 

^2pki\vk) 

k=l 



(36) 



where p k = \p [k mod K) . 

Second, because the blocks D in Eq. (|A-9[) in the Ap- 
pendix are the same, we note that all diagonal elements 
Da from re must be zero, otherwise only t = will satisfy 
7(g) — tKA © kb > 0. This means that the only states, 
which can appear in the expansion (|36[) are of the form 
\rik) = |fcfc), since the \a k ) and \b k ) have to be eigenstates 
of the operators Di = \i)(i\ (see the Appendix A). ■ 

Having proved the optimality of the expansion of ka © 
res in Eq. (|30| we can now provide the general formula 
for the function £ (\ip)) for pure states in the Schmidt 
decomposition. The value of the function V(\ip}) is given 
by the solution of the following max-min problem 

a = max min ~ V 7 ^) ; 1 < i < j < d, (37) 
V i<3 Pi + Pj 

where the first max is taken over all possible probability 
distributions V = {pi,P2-, ■■■}• A solution of this problem 
for the case d = 3 and d — 4 is given in Appendix and 
we can summarize: 

Proposition 11 (£ for pure states), (a) If \ip) = 
Si=i V^iN) i s a pure two-qutrit state, then 



£(\^)) = 2^X~X~ + 2^X~XT - X io , 



(38) 



where lojioi^o are pairwise different and jo,k are such 
that (y/Xfo- \fhT Q ) 2 > ( V^i - V^) 2 f or al1 3, k - 
(b) If\ip) = 2i=i V^i|w) * s a pure state in a AxA-system, 
then 



£{\ip}) = max{ei,e 2 ,e 3 }, 



(39) 



where 



ei = 2\J AiA 2 + 2\J A3A4, t2 — 2yJ A1A3 + 2\J A2A4, 

e 3 = 2 v /AiX+2 v / A^A 3 ". (40) 

Note that in both cases we have for a maximally en- 
tangled state £(il>) = 1- 
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Schmidt-correlated states 



which leads to the following max-min problem for V(gsc) 



To conclude the section we consider a family of mixed 
states, for which the introduced function £{q) can be 
also computed exactly. These states are called Schmidt- 
correlated (SC) states in the literature [l(|. By defini- 
tion, SC states are a mixture of states that share the 
same Schmidt basis 



N 



Qsc = with 



u=i 



(41) 
(42) 



i=i 



SC states can be written in computational basis directly 
as 



T/ / v . Qu + Qjj - %Qij 
V(qsc) = max mm — - 

V i<j Pi + Pj 



max mm ■ 



Pi + Pj 



(47) 



This problem can be effectively solved numerically or 
with the methods of the Appendix B and its solution 
gives the exact value of the function £(gsc)- For two 
qubits one finds 



(48) 



as a nice analytical expression. 



Qsc 



with Q k 



,(«) 



(43) 



As in the case of pure states, we find for the SC states 
the optimal decomposition of ka © ■ 

Lemma 12 (Optimality for SC states). In the case of 
SC states the optimal decomposition of ka © kb for the 
estimation of the parameter £{g) can always be written 
in the form of Eq. i3(A) : 



V. THE ENTANGLEMENT PARAMETER £{g) 
AS A LOWER BOUND ON THE CONCURRENCE 

In this section we demonstrate that the function £(g) 
can be used to estimate the amount of entanglement in a 
quantum state. More specifically, we show how it delivers 
a lower bound on the concurrence, which is a well known 
measure of bipartite entanglement. For bipartite pure 
states in a d x <i-system the concurrence is defined as 
EE El: 



opt 



© K 



opt 



Pil(\H))- 



(44) 



cm) 



i 



l-Tr(^). 



(49) 



Proof: There were two essential ingredients in the 
proof of the Lemma [lOj First, we used the fact that 
the CM of a state, written in Schmidt decomposition, 
is invariant under interchange of parties. Obviously the 
same invariance does also hold for SC states. Second, we 
used the fact, that all blocks D of the CM are the same. 
Using the formulae of the Appendix A one easily verifies 



that Dg C 



u sc 



u sc- 



For these states the problem of calculating the function 
£{q) reduces to the max-min problem in Eq. ([37)1 . This is 
due to the fact that diagonal elements of the covariancc 
matrix have a pretty simple form for qsc- Indeed, using 
the formulae from the Appendix A we calculate directly: 



(Df c B/C ) 

A/B 



X SC 

yC 
A 5C 



Y, 



A/1 



SC 

I sc 



Qnkj - QuQjj, 1 < i < d, 
3 = idiag{^ + Qkk}, 1 < i < k < d, 



diag{ft fe }, 1 < i < k < d. 



(45) 



The 2x2 blocks in Eq. d3TJ) will then take the form 



n 2x2 



Qii + Qjj - t{pi+Pj) 

±2Qij 



±2g t j 
Qii + Qjj - t(pi +Pj) 



(46) 



In this definition, we introduced already a prcfactor 
which guarantees that < C < 1, this will turn out 
to be useful for our purposes. 

The concurrence is then extended to mixed states by 
the convex-roof construction 



C(Q) 



Pi:\"4>i) 



PiC(|^)), 



(50) 



where the minimization is taken over all possible decom- 
positions of the state g = EiPilV'i)(V'j|- Of course, this 
minimization is quite difficult to perform, and only for 
two-qubits a complete solution is known [18| . Therefore, 
it is desirable to have at least some lower bounds on the 
concurrence. 

The idea of obtaining lower bounds on C from £ is as 
follows: Let us assume that one can prove a lower bound 
like 



C(M) >a£(m+P 



(51) 



for pure states only with some constants a, (3 and a > 0. 
Then, since £ is convex, the right hand side of Eq. (fBTj) 
is convex, too. By definition, the convex roof is the 
largest convex function which coincides with C on the 
pure states. Consequently, C(g) > a£(g) + (3 holds for 
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all mixed states, too. This trick has already been em- 
ployed in several works to obtain lower bounds on entan- 
glement measures [3, However, as the CMC detects 
many bound entangled states where other criteria fail 
plj . our results will deliver entanglement estimates for 
states, where the other methods fail. 



A. Two qubits 

Using the Schmidt decomposition, one can express the 
concurrence for pure states in terms of Schmidt coeffi- 
cients as 



cm) = 



2d 



i<j 



(52) 



Comparing Eq. ([52]) and Eq. ([29]) from the Section HVl we 
see that the concurrence and the function E(\tp)) coincide 
for on two-qubit pure states 



Consequently, C{g) > £(g) holds for any mixed state. 
Note, however, that for the special case of two qubits 
one can calculate the concurrence also directly for mixed 
states HI. 



B. Two qutrits 



Using the solution of the problem (|37|) it is possible to 
derive a lower bound on concurrence for pure states of 
two d-level systems. Before we proceed, note that j3] 



cm = v d^W^> Aj - dzrE v^7. ( 54 ) 

y i<j i<j 

This follows from the fact that 

E A ^ = duhi) E E^ + A ^<) ( 55 ) 

i<j ^ ' i<j k<l 



d(d-l) 



i<j k<l ' i<j 

For two qutrits £(\ip)) is given by Eq. ([55]) . We have that 
2y/XiXj + 2\J XiXk — Xi = 2y/XiXj + 2\f~XiXk + 2^/XjXk 



— 2yJ XjXk — 1 + Xj + Afc 
<2C(|^)) + ( V ^7- V ^) 2 -1 
< 2C(|^». 

Hence we have for mixed two-qutrit states 



(56) 



(57) 



Using the results from Proposition [6] we have, for in- 
stance, 



C(Q) > 



lf Tr(gj) + Tr(g|)-2 



4 I 



+ ^[Tr(^)-Tr(^)P + ||C||? r }, (58) 



which is an easily computable lower bound that delivers 
non-trivial estimates for many weakly entangled states. 



C. 4x4 systems 

In this case £ (\ip)) is given by Eq. (|39|) . Wc can directly 
estimate: 



£{\ij})) =max{ei,e 2 ,e 3 } < 20^+2 ^X~ 4 (59) 
+ 2^AiA 2 + 2VA3A4 + 2VA1A4 + 20^ 



(53) and hence for arbitrary mixed states 



C(g) > f(g). 



(60) 



D. Examples 

Let us discuss the strength of these lower bounds by 
considering some examples. Let us first consider Bell- 
diagonal two-qubit states. For them, the reduced states 
qa and qb arc maximally mixed, and then Proposition [6] 
delivers the bound C(g) > Tr(|C|) - 1/2. On the other 
hand, it is known that for Bell diagonal states the con- 
currence is given by C(g) = 2A max — 1, where A max is the 
maximal eigenvalue, i.e., the maximal overlap with some 
Bell state [l7[- Noting that A max = [1 + 2Tr(|C|)]/4 (this 
can be easily seen if the closest Bell state is the singlet 
state and we take appropriately normalized Pauli ma- 
trices as obscrvables in the definition of the matrix C), 
one finds that our lower bound is tight for Bell diagonal 
states. 

For general two-qubit states, the lower bound cannot 
be tight, as they are entangled two-qubit states, which 
are not detected by the CMC. On the other hand, any full 
rank two qubit state can be brought to a Bell-diagonal 
state by filtering operations. Since it is known how the 
concurrence changes under filtering operations (2p| . one 
could use the filtering and our lower bound to determine 
the concurrence for arbitrary two-qubit states. 

For two qutrits, our bound is not tight for states like 
\i>) = (|00) + |11) + |22»/V3 or = (|00) + |22))/V2, 
however, for the latter the reason lies in the fact that 
the bound ([54]) is not tight. On the other hand, the pre- 
sented method delivers nontrivial lower bounds for many 
bound entangled states (such as the the family of chess- 
board states), as many states of this type are detected 
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by the CMC [Tl|, but not by the PPT or CCNR crite- 
rion (which means that the methods from Ref. [3| must 
fail). Similarly, our methods can be used to estimate the 
entanglement of bound entangled states for 4 x 4-systems. 



VI. CONCLUSION 

In conclusion, we have introduced an entanglement pa- 
rameter £ that quantifies the violation of the covariancc 
matrix criterion. We have shown that this parameter is 
convex and invariant under local rotations, but it can 
increase on average under local operations and classical 
communication. Most importantly the parameter £ can 
be used to deliver lower bounds on the concurrence. 

For future work, it would be interesting to connect £ to 
other entanglement measures, such as the entanglement 
of formation Even more interesting, would be an 
extension of the covariance matrix criterion to the mul- 
tipartite case and a definition of a similar entanglement 
parameter there. This could help to quantify entangle- 
ment in multipartite systems, where much less is known 
compared to bipartite systems. 

We thank Jens Eisert, Bastian Jungnitsch, Matthias 
Kleinmann, and Sonkc Nickamp for discussions. Espe- 
cially we thank Philipp Hyllus for discussions and com- 
ments on the manuscript. This work has been supported 
by the FWF (START Prize) and the EU (OLAQUI, 
QICS, SCALA). 



APPENDIX A. 

Here we calculate symmetric block CM of a pure bipar- 
tite state, which is written in the Schmidt decomposition 
|V') = yi f -v/ATIm) ® \ib)- Consider d A = d B = d. As 
it is proven in [1J| we can choose the basis in the oper- 
ator spaces B(Ha) and B(Hb) arbitrarily for applying 
the CMC. In this case it is convenient to choose the local 
orthogonal observables 



(A-l) 



The mean values for the state 



X itj = -j= + \mi l<i<j<d, (A-2) 



A = i = l,...,d, 

1 

Y k ,i = -j=(\k)(l\-\l)(k\), l<k<l<d 1 (A-3) 

which satisfy following anticommutation relations: 
{D h Dj} = 5 t] + , {A, Xij} = Xu, 

{ &i i } — Y\ j , { Xi j , Yij } , 

{Xij, X^} = Di + Dj, Yij} = Di + Dj. 

(A-4) 

Note that this is not the complete set of relations, how- 
ever other relations will not give any contribution to the 
CM and hence we leave them out here. 



(*£®i> = <i®x£> = 0#®i 

(D A ®1) = (1® D B ) — A l 
{{D?,Df}® 1) = (Aj + 



are given by 
. = (1 



3tf> = 



The blocks A, B and C of the 7 (IV')) can be therefore 
written as 3 x 3 block matrices. Because of the relations 
(|A-4|) and (|A-5|) a lot of terms in these blocks will be 
equal to zero an we have the structure 

/ D A/B/C o \ 

A,B,C=\ x A ' B l c , (A-6) 
V yA/Bic J 

Since the off-diagonal terms can be calculated straight- 
forward 



(Di 



(Dt)(D] 



XiSij 



AjAj, 



>Y£) = (X±®Y r B ,)=0, 



(X pq (g) X rs ) — \J XpXqSprSqg, 
(Y p ^ (g) Y B S ) = —y/XpXgSprSqs, 

we can write the blocks in (|A-6j) as follows 



(A-7) 



D = 


D A/B/C 


ij 


X 


= X A ' B 


Y 


= Y A I B 




X c 




yC 



— XjSij — XiXj, 

= -diag{A 4 + A fe }, 1 < i < k < d, 
= -diag{Ai + A fe }, 1 < i < k < d, 



= diag{ y/XpXq}, 1 < p < q < d, 

= diag{- v /V^},l <P< q<d. (A-8) 

Finally, we arrive at the general form of the CM for a 
pure state as a function its Schmidt coefficients: 



(D D \ 

X X c 

Y Y c 

DO D 

X c X 

Vo or c o o y J 



(A-9) 



with the blocks given in Eq. (|A-8|) . 

APPENDIX B. 

In this Appendix we discuss the possible ways of solv- 
ing the max-min problem: 



t = max min ■ 



-, 1 < i < j < d. (B-l) 



T i<3 Pi + Pj 

We consider the cases d = 3 and d = 4. We define 

b l3 = (Va7- \^)\ (B-2) 



Pi + Pj 



(B-3) 
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1 2 3 1 2 3 4 




FIG. B-l: (a) In the case d = 3 there are only three elements 
a.ij. These elements are written in the form of tableaux in 
order to visualize the problem considered, (b) Tableaux of 
the elements in the max-min problem for d = 4. 

For d = 3 there are only three different a's that can be 
arranged in a tableaux as in Fig. IB-lf a) . 

The properties of the solution can be summarized as 
follows: 

Lemma B-l. (a) Consider the optimization problem in 
Eq. &B-1\) for d=3 with the only assumption that bij > 0. 
Let jo , ko be such that 

b joko = max{6j7c} (B-4) 

Then the optimal solution oP is given by 

a = min {a 1 , a 11 } , (B-5) 

where 

a 1 = ^(612+613 + 623) (B-6) 
a " = Kdo +Km (B-7) 
with j ^ i ^ ko- 

(b) For the same problem, if the bij are given via 
Eq. W-2]) as functions of Schmidt coefficients and ful- 
fill therefore further restrictions, the optimum is always 
given by 

a° = a 11 = 1 + A l0 - 2 V /A~A~ - 2^/X la X ko (B-8) 

Then we also have that a 11 = min.yfe{l + A.; — — 
2\/\i\k} where the i,j, k are pairwise different. 

Proof: (a) Let us first assume only that 6y > 0. In the 
max-min problem (|B-ip the maximization is taken over 
all possible probability distributions. It is convenient to 
distinguish two cases: 

Case 1: The optimal probability distribution does 
not have any zero elements. Assume Vo is the opti- 
mal distribution and p\ ^ Vz. We often drop the 
index in the following for simplicity. We show that 
this optimal distribution necessarily has to be such that 
d\2 = o*i3 = «23 = or, otherwise the optimality is vio- 
lated. Indeed, assume that this is not case. Then, with- 
out loss of generality, we can write an < «i3 < a23, 



where one of the inequalities must be strict. Now con- 
sider some distribution V' such that 

p'l =Pi- 2e, p' 2 = P2 + e, p' z = P3 + e, (B-9) 

with some e > 0. The coefficients will change and 
become according to the new distribution V' 

a'l2 > "12, "13 > "13, Ck'23 < Qf23- (B-10) 

Since the parameter e can be chosen arbitrarily small 
the number a' 12 will be still the minimal one, i.e. a' 12 — 
mini<j ajj- But a[ 2 > «i2- Consequently the distribu- 
tion V' gives a bigger minimum of the set {o^ } than the 
distribution Vo, which contradicts the assumption that 
Vo is optimal. Hence we conclude that a\ 2 = 023 must 
hold, which implies ai2 = "13 = Q23- 

Having established that if Vo is optimal and contains 
no zero elements, then a\<x — a\s = "23 = a holds, we 
can calculate a explicitly. We have 



Pi +P2 Pl+ P3 P2 + P3 

Multiplying by the denominators summing up these 
equations gives 

2a°( Pl +p 2 + pa) = 612 + 613 + 6 23 . (B-12) 
Because p\ + p 2 + P3 — 1 we arrive at 

a / = a° = i(6i2+6i3 + 6 23 ). (B-13) 

Case 2: The optimal probability distribution Vo has 
at least one zero element. This means that one ctij = 6y 
independently of the two free parameters of the proba- 
bility distribution (since pi +pj = 1). We can distinguish 
three cases, and assume for defmiteness 6 i2 < 613 < 623- 

(i) If ai2 = 612 (that is, P3 = 0), then clearly > bij 
for i,j = 1,3 and i,j = 2,3. Then we have min{ay } = 
o;i2. But then decreasing one of the pi or p 2 and increas- 
ing consequently P3 will lead to an increasing of 012 and 
a better solution which belongs to case 1. So a solution 
with ai2 = 612 can never be optimal. 

(ii) If q;i3 = 613 the optimal probability distribution 
has to be such that a\s < ct\2 and ai3 < 023. But as in 
the case (i) one can directly see that this leads to case 1 
and can never be optimal. 

(iii) Finally, consider the case a.23 = 623. Then, one can 
see as in case 1 one can achieve a\2 = a\3 without giving 
up optimality. More precisely, the optimal probability 
distribution has to fulfill this from the beginning (if ai2 
and Q13 are the minima) or it can be achieved (if 023 is 
the minimum). 

This leads as in case 1 to the conclusion, that we have 

ai2 = — = — => ai2 = 612 + 613, (B-14) 

P2 P3 
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and consequently 
bi2 



a 



II 



"12 



13 = 1 + Ai - 2^M\ 2 - 2VA1A3. 

_ (B-15) 

However, it is not yet clear what the min{o!y} is. Two 
cases can be distinguished: 

(iiia) If ai2 > OL23, — 623 one would take minjay} = 
"23 = ^23 j but then, one can improve it further as in 
the cases (i) and (ii) by going to the case I and taking 
finally a 1 from Eq. (|B7L3"]) . Note that a 1 = {a 11 + b 2 z)/2. 
Therefore, if a 11 = ai2 > "23 = &23 one has also that 
a 1 < a 11 , so effectively one takes min{a 7 , a 11 }. 

(iiib) If Q12 < CK23 = 623 we take minjajj} — a 11 and 
going to case 1 does not help. But in this case, we have 
a 1 > a 11 , so effectively one takes again minja 7 , a 11 }. 

Finally, let us discuss shortly the meaning of the choice 
jo and fco in Eq. (|B-4|) as one may consider also in 
Eq. (|B-7|) with other indices. However, one can directly 
compute that aj^ < a 1 is equivalent to 6 ?J + bn~ < bjk 
and this can only be true, if j and k are chosen as in 
Eq. (|B-4[) . In other words, the a 1 ^ for other indices than 
jo,ko can never contribute and one could alternatively 
write that ao = minja/, a[ 2 , a 2 ^}. 

(b) Let us now assume that the bij stem from Schmidt 
coefficients as in Eq. (|B-2[) . We know from the previous 
discussion that we have to take a 11 iff bi j +bi k a < bj k a ■ 
In terms of the Schmidt coefficients, this implies that 



Akor > (va 



A J0 r + G 



- V Afc ) ■ 
(B-16) 

This, however, is true for any triple of positive real num- 
bers vXJ, if Jo and ko are chosen as in Eq. (|B-4j) . Then, 
its also clear that the a 11 chosen is minimal among all 
thebij+bik- ■ 

Further, we discuss the case d = 4. The elements a,-j 
are again embedded in a tableaux as in Fig. IB-lf b). We 
begin with studying of properties of the optimal proba- 
bility distribution Vq. Suppose as in the case d = 3 that 
ctij correspond to the optimal probability distribution Vo 
and that a\ 2 



miny{a? }. We can formulate: 



Lemma B-2. The solution of the max-min problem |5 
for d — 4 is given by 



a 



°-min{aVV m 



}• 



(B-17) 



where 



a 1 = l-2^/xJ^-2^/X 3 ~)^, 
a 11 = l-2y/\iM- 2v^Al, 



jii 



1 -2^/x^x i -2y/x 2 ~x 3 . 



(B-18) 



Proof: The proof proceeds in several steps. 
Step 1. Let us first consider optimal probability distri- 
butions Vo = {pi,P2,P3,P4} where all pi are nonzero. In 
this case we show that for a° = min{ajj} at least one of 



the three equations must hold 



a ' 



"12 = "34, 
"13 = "24 
"14 = "23- 



(B-19) 



The idea of the proof is similar to the proof of Lemma 
IB-It we consider small perturbations of the optimal prob- 
ability distribution Vo that increase the minimal clement 
a and therefore destroy the optimality if some addi- 
tional constraints are not fulfilled. As we will see, these 
constraints will give us the conditions Eq. (|B-19p . 

Let us assume for definiteness that the optimal a is 
given by ai 2 . We can consider the following four trans- 
formations of the pi : 



T\ : 


p'i 


= Pl 


- 3e, 


P'i 


= P% 


+ e 


for 




T 2 ■ 


P-2 


= P2 


- 3e, 


P'i 


= Pi 


+ e 


for 


i £ 2, 


T 3 : 


P's 


= P3 


+ 3e, 


P'i 


= Pi 


— e 


for 


i ^ 3, 


T 4 : 


P'4 


= P4 


+ 3e, 


P'i 


= Pi 


— e 


for 


i ± 4, (B 



where e can be chosen arbitrarily small. All the trans- 
formations increase a\2, but all have to keep the op- 
timality of the probability distribution, so that mini- 
mal a given by V' cannot be larger than u\2- From 
transformation T\ it follows that Vo is optimal if and 
only if CJ12 = min{o!23, "24, "34}, as these entries de- 
crease under the transformation. Similarly, it follows 
from T2 that a.\ 2 = min{ai3, an, "34}, and from T3 
that ai2 = min{ai3, "23, "34}, and finally from T4 that 
a\2 = min{ai4, Q24) "34}- Given this finite number of 
possibilities, one can directly check that cither ai2 = "34 
or CK12 = "13 = "24 or a>i2 = "14 = "23 must hold for 
optimal probability distribution Vq which proves the first 
claim. 

From these conditions we see that there are the three 
candidates for the optimal a : 



n° 


= "12 


= "34, 






= a ! = l 


a 


= "13 


= "24, 






= a n = 


n° 


= "14 


= "23, 






= a 111 = 



1>V2 



1 



;i3 



2^X1X2 — 2\J A3A4, 
6 24 = 1-2^*1^-2^^1, 



>i4 + &23 = l-2 v /A 1 A 4 "-2v / A 2 A^. 

(B-21) 

Step 2. At this point, we have identified three candi- 
dates for the a , but is is not clear yet, which one should 
be taken. 

We will show now, however, that only the minimum of 
these can give a valid solution. For that, assume that one 
has a probability distribution V\ which has the optimal 
a°(Vi) = a 1 . Then a® 4 = a\ 2 = min^ay} and hence 

"12 < "13 =>■ hzipi +P3) < bia(pi +P2), 

"34 < "23 =>■ b 3A (p 2 +P4) < b 2i {p3 + Pi)- (B-22) 
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Consequently, 612 + 634 < 613 + 624 and hence a 1 < a 11 . 
Similarly, it follows that o 7 < a 111 . So if one finds a 
solution, then it has to be the minimum of all a k . 

This also shows that if there is a second solution V 2 
with a°(P 2 ) = a 11 , then a°(Pi) = a°(P 2 ) must hold, 
since a 1 < a 11 and a 11 < a 1 . Note also that the ar- 
guments leading to this did not require the assumption 
that the probability distributions have nonzero elements. 

Summarizing Step 1 and Step 2, we can state that if 
there is a optimal probability distribution with non-zero 
elements, then the solution is given by 

a =mm{a I ,a II ,a 111 }. (B-23) 

Step 3. Now we have to consider the cases where the 
optimal probability distribution has some zero elements. 
Let us first consider the case that there is exactly one 
zero element. 

There exist two possibilities. The first one arises, when 
the minimum is given by ai2 and p^ = 0. Then, the 
transformations T\ , T 2 and T4 in Eq. (|B-20|) can still be 
applied, but we have to modify T3, since there are no 
negative probabilities 

^3 : P3 = z>3 + 2e, p'i =Pi-e for i = 1,2, p' 4 = p 4 = 0. 

(B-24) 

This transformation leads exactly to the same condition 
as T3 above ai% = min{ai3, 023, 0(34}. Therefore, the 
same conclusion as in Step 1 can be drawn. Similarly, by 
considering T4, one can show that if p 3 = 0, the conclu- 
sion from Step 1 still holds. 

The second possibility arises, if the minimum is again 
given by ai2, but this time p\ = 0. Then, only T 2 in 
Eq. (|B-20j) can be applied. We define the modified trans- 
formations: 

T3 :p 3 = P3 + 2e, p'i =Pi-e for i = 2, 4, p\ = pi ; 

T 4 :p4=P4 + 2e, p\ = Pi - e for i = 2, 3, p[ = pv, 

(B-25) 

Then, repeating the argumentation from Step 1, one ar- 
rives at the same conclusion, apart from the special case: 
«12 = ai3 = «i4 < a>ki for k,l £ {2,3,4} holds. 



In this special case, we have that a = (bi 2 + 613 + 644) 
and consequently = &ifc/(&i2 + fri3 + &i4) for k = 2,3,4. 
Since a 23 = b 23 /(p 2 +p 3 ) > a = (612+613+^14) it follows 
that 623 > bi 2 + 613. Generally we have bki > &i& + b\i, 
for k,l £ {2,3,4}. 

Due to the definition of the by, it means that the 
Schmidt coefficients have to fulfill 

{^k-V\if > (v / Ai"-Y / A^) 2 + (v / Ai"-Y / A7) 2 (B-26) 

for k, I £ {2,3,4}. Since the V^fe are positive real num- 
bers, this can only hold if y/Xi is inside the interval 
[V^k ; \fXi\ ■ As there are three intervals, and two of them 
intersect in only one point, we must have that y/X[ = V^7 
for some i £ {2, 3, 4}, which implies that the correspond- 
ing b\i = and an = 0. Since a\ 2 = ais = an all of 
them must be zero and hence a = and all b^ = for 
any k £ {2, 3, 4}. Physically, this means that all Schmidt 
coefficients are the same and the state is a maximally en- 
tangled one. But then also a 7 = a 11 = a 111 = 0, so this 
special case does not deliver a novel solution. 

Step 4- Let us now consider the case, where two or 
more pi equal zero. 

Let us first assume that exactly two pi are zero, namely 
P2 = P3 = 0. Then an = 614 and 023 = 00 are indepen- 
dent of the probability distribution. However, if we make 
the transformation 

1 : pi = Pi - £ i £ {1,4}, p' k = Pk + e k £ {2, 3}, 

(B-27) 

the minimal value a does not decrease (as all 
ot\2, «13, Q!24, a34 remain constant and ai4 increases). 
Therefore we arrive at a solution, where none of the pi is 
zero and which is as good as a solution with p 2 = p% = 0. 
Thus we conclude that solutions given by distributions 
with two zero elements are contained in solutions char- 
acterized in Step 1. 

Finally, we have to discuss the case that three pi equal 
zero and consequently the remaining one equals one. 
This can be excluded with a similar transformation as 
in Eq. (|B-27[) and we leave the details as an exercise to 
the reader. ■ 
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